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ABSTRACT 

Radial velocity surveys find Jupiter mass planets with semi-major axes a less than 0.1 AU around 
~ 1% of solar-type sta rs; counting planets with a as large as 5 AU, the fi'action of stars having 
planets reaches ~ 10% (|Marcv et al.ll200"5t iButler et al.ll2006f) . An examination of the distribution of 
semi-major axes shows that there is a clear excess of planets with orbital periods around 3 or 4 days, 
corresponding to a « 0.03 AU, with a sharp cutoff at shorter periods (see Figure 1). It is believed 
that Jupiter mass planets form at large distances from their parent stars; some fraction then migrate 
in to produce the short period objects. We argue that a significant fraction of the 'hot Jupiters' 
(a < O.IAU) may arise in binary star systems in which the orbit of the binary is highly inclined to 
the orbit of the planet. Mutual torques between the two orbits drive down the minimum separation 
or periapse rp between the planet and its host star (the Kozai mechanism). This periapse collapse 
is halted when tidal friction on the planet circularizes the orbit faster than Kozai torque can excite 
it. The same friction then circularizes the planet orbit, producing hot Jupiters with the peak of the 
semimajor axis distribution lying around 3 days. For the observed distributions of binary separation, 
eccentricity and mass ratio, roughly 2.5% of planets with initial semimajor axis « 5 AU will migrate 
to within 0.1 AU of their parent star. Kozai migration could account for 10% or more of the observed 
hot Jupiters. 

Subject headings: binaries:general;planetary systems;celestial mechanics 



1. INTRODUCTION TO KOZAI MIGRATION 

Statistics f rom r adial velocit y pla net searches 
(|Marcv et all 120051 iButler et all I2006D show that 
the occurrence rate of giant planets within 0.1 AU 
( "hot- Jupiters" ) is ~ 1%; extrapolating to 20 AU the 
occurrence is 12%. There is a clear "pile- up" of planets 
with orbital periods near 3 days (Fig. [T|). Transit 
observations yield a similar fraction of hot Juptiers 
(|Gould et al.ll2006HFressin et al.l[200l . What migration 
mechanisms can produce the observed feature in semi- 
major axis distributions represented by hot Jupiters? In 
this article we focus on the mechanism known as Kozai 
migration. 

Consider a planet circling a star that is a member of 
a binary system. The mutual torques between the bi- 
nary and planetary orbits transfer angular momentum 
between the two while leaving the orbital energies nearly 
unchanged. For mutual inclinations / ^ 40° a resonance 
between the precession rate of the planet's nodal and 
apsidal lines greatly enhances the effectiveness of this 
exchange of angular momentum, producing large oscilla- 
tions in_the planet's angular momentum (Kozai cycles, 
lKozai|[r962l ). The planet eccentricity (cp) and periapse 

i rp = ap(l — Bp)) osc illate with a characteristic timescale 
Holman et~aLlll997D 
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where and rric are the masses of the central and com- 
panion stars, while Pc and Pp are the periods of the 
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binary and planetary orbits, respectively. The b inary 
eccentricity is denoted by Cc. iHolman et al.l (|1997f ) and 
iTakeda fc Rasiol (|200l . among others, have studied the 
role of these Kozai cycles in producing the eccentricities 
observed in known exo-planets. 

For sufficiently large /, rp can reach very small values, 
allowing tidal dissipation to erode the orbit of the planet. 
lEggleton fc Kiseleva-EggletonI ()2001[ ) were the first to 
propose that Kozai cycles, in combination with tidal fric- 
tion, can shrink the orbit of a inner binary in a hierar- 
chical trip le system, lead i ng to the formation of contact 
binaries. IWu fc Murravl (|2003l . hereafter WM03) have 
studied Kozai migration in application to exo-planets 
and found it to be the only plausible explanation for the 
migration of the planetary object IID80606b. 

In the absence of any other modification of the gravi- 
tational potential, the minimum may fall below the 
Roche radius (r/j) and the planet may be destroyed. 
However, there are a number of competing torques that 
can limit the amount of angular momentum that the 
Kozai torque can extract from the orbit of the planet, 
including general relativistic (GR) corrections to Newto- 
nian gravity, and torques associated with the extended 
mass distribution of both the primary star and the 
planet. The latter includes rotationally induced plan- 
etary oblateness, the tidal bulge raised by the star on 
the planet, the misalignment of this bulge produced by 
friction, and the stellar counterparts of all these. These 
torques can halt the Kozai-induced collapse in Tp and 
promote planetary survival. 

Which torque becomes competitive with the Kozai 
torque depends on the system; for systems with very 
large binary semi-major axis (oc) and therefore very weak 
Kozai torque, the GR precession can halt the reduction 
of rp before tides become important. However, for tighter 
or more inclined binaries, tidal friction sets the minimum 
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Tp. Since the tidal torques depend strongly on rp, binary 
systems with a wide range of ad ap will be stalled at es- 
sentially the same r^, leading to a pile-up of hot Jupiters 
at Up ~ 2rp when the planet orbits are later circularized. 

2. NUMERICAL EXPERIMENTS 

We quantify the effect of Kozai migration by consid- 
ering an ensemble of binary systems following that in 
iTakeda fc Rasiol (|2005| ) . These binaries are initially com- 
prised of a solar-mass host star, a jupiter-mass planet 
{nip = Mj) orbiting at 5 AU with an eccentricity of 
0.05, and a binary companion of mass 0.23Mq - this is 
the peak of the obser ved mass ratio distribution in the 
solar neighbourhood (|Duauennov fc Mavoij|199il ). The 
distribution in binary separation (P(ac)) is assumed flat 
in logarithmic Oc (cc ranging from 20 to 20, 000 AU). We 
set P{eb) — 2ef,, a thermal distribution often adopted in 
binary population synthesis. This latter choice hardly 
affects the results. The last ansatz, our most sensitive 
yet most uncertain assumption, takes / to be isotropi- 

cally distributed; B ased o n studies of stellar spin and 

binary orbits (jHald |1994( ). this seems reasonable for 
ttc > 40 AU, but may be less appropriate for tighter 
binaries; polarimetry studies of protostellar disks sug- 
gest that the circumstellar disk and the binary pl ane are 
corre l ated for up to a few hundred AU ( Jense n~et al.l 
12004 iMonin et alt I2006D . However, polarimetry esti- 
mates only the projected angle between the two planes, 
and is strongly plagued by interstellar polarization. The 
results should be taken with caution at present. 

We produce an ensemble of 100, 000 systems. Out of 
these we select systems that can potentially perturb the 
planet to Vp <^ 0.1 AU. To reach this distance, a planet 
starting at semi-major axis Op (with a small eccentric- 
ity) will have to attain Cmax > 1 — O.l/cp. Ignoring 
tidal dissipation,'* the Kozai integral (the planet's or- 
bital angular momentum in the normal of the binary 

plane) Hk = 



'(1 — ep)cos/ = constant. Taking a 

minimum / » 40deg during the Kozai cycles (see, e.g., 
iHolman et al.| [l997'l. this yields a minimum initial incli- 
nation required for producing hot Jupiters: 1^81 deg. 
This value is independent of the binary separation or 
mass. The fraction of isotropically inclined systems that 
have such a misalignment is ^ 15%. 

We then weed out planets that are likely dynamically 
unstable according to the following fitting formula 



> 0.330 

an 



Allen + 0.069e, 



(2) 



This expression is obtained by integrating the orbits 
of our initial system for 10^ binary orbits, taking / = 
90 deg. This non-coplanar stability limit is 15% to 30% 
m ore restrictive than t he coplanar stability limit found 
bv IHolman fc WiegertI {1999). It is used here as a rough 
proxy for systems that either eject their planets quickly 
after formation, or are unable to form planets due to the 
strong tidal influence of the companion. This proceedure 
eliminates many systems with Oc < 100 AU; we are left 
with ^ 10% of the original ensemble that could poten- 
tially reach < 0.1 AU, if they arc not stalled by other 
torques at larger distances. 

^ Tidal dissipation increases the Kozai integral and slightly raises 
the minimum requirement on / (WM03). 
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Fig. 1. — The histogram of the planet semi-major axis (loga- 
rithmic) distribution. The thick solid curve is the observed radial 
velocity planet distribution. The planets detected by the transit 
technique are added on top (dotted line) assuming that the de- 
tection efficiencies are the same between the two techniques. The 
shaded area shows the simulation result, with the vertical axis read 
at the right. The peak at ~ 3 day orbital period corresponds to 
planets that are Kozai migrated and later circularized. The posi- 
tion and width of this peak depends on a number of parameters 
(see Eq. [6]). In particular, if planet radius is a decreasing function 
of planet age, the width of the peak shrinks (see Fig. [2]l. The nar- 
row peak at ^ 5 AU corresponds to planets that are unmigrated, 
remaining at their initial ap. 

These remaining systems are integrated using secu- 
lar equations obtained by averaging over the orbital 
motions of both the planet and the binary companion 
(|Eggleton fc Kiseleva-Eggletonll200if) . These equations 
include the effects of Kozai perturbation, tidal dissipa- 
tion, GR precession, and tidal and rotational bulge pre- 
cessions.^ We use a Runge-Kutta integrator with an 
adaptive step size set to keep the integration error be- 
low a preset limit. We follow the procedure described in 
WM03, which also lists values for the various parameters 
involved. In particular, we choose the initial stellar spin 
direction to be aligned with the initial orbit normal for 
the planet. 

The integration is stopped after 5 Gyrs have passed, 
or when 5 million timesteps are exhausted, or when 
rp < 2Rq . The last condition roughly corresponds to the 
planet overflowing its Roche lobe; however, none of the 
planet in our simulation reached this state. ^ The limit 
on the number of integration timesteps is usually reached 
if Kozai oscillations have been effectively halted by rapid 
tidal or other precessions; in that case the subsequent 
dynamical evolution of the planet simply reduces Cp. We 
then use a simplified code, including only the effects of 
tidal dissipation on the planet orbit and planet/stellar 

^ In this study, we rely exclusively on these secular equations. 
The actual dynamics may deviate due to short term noises and 
mean-motion perturbations and should be studied with N-body 
integration codes. 

® This is due to the strong dependence of the tidal timescale on 
rp] tidal distortions act as a barrier, maintaining rp >rii. 
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Fig. 2. — Similar to Fig. [T] except where we have taken the 
planet radius to shrink as Rp = Rj[l + exp{— 4/3 X lO^j/rs)]. The 
3-day feature narrows significantly. 

spins, to finish integrating to 5 Gyrs. 

We find that about 2.5% of our ensemble eventuaUy 
migrate inward of 0.1 AU. The distribution of final semi- 
major axes is concentrated between 0.02 AU and 0.05 AU 
with a peak at 0.03 AU. Our hot Jupiters exhibit a pile- 
up at ~ 3 day periods similar to the observed population 
(Fig. HD. 

Given the same initial /, tighter binaries produce a 
closer-in hot Jupiter in a shorter amount of time. Many 
of the hot Jupiters are tidally ensnared on their first close 
approach to the host star (Fig. [3]), with the Kozai period 
between 10^ to 10® yrs. Tidal circularization of these 
orbits then takes upward of 10^ years. 

The 3-day feature in the computed Op distribution ap- 
pears wider than the observed distribution. However, 
as just noted, closer-in planets are migrated in ear- 
lier, so they still have larger radii and larger stalling 
peraipses. Experimenting with the following time evo- 
lution of planet radius. 



exp 



t 



'^shrink 



(3) 



with Tshrink taken to be 3 x 10^ yrs, we find that the 3-day 
bump narrows significantly (Fig. [2]). 

3. DISCUSSION 

3.1. Stalling Radius and the 3-day Pile-up 

The periapse of a Kozai-migrating planet is stalled 
at a distance where the eccentricity forcing due to 
the binary companion is counteracted by the eccentric- 
ity damping by tidal d issipation. Kozai forcing yields 
(jEggleton fc Kiseleva^ gglcton 2001.) 

1 dep _ ^^-^ f^cnl 1 

ep dt ^ "P mp + m, -I- /^n^^l ~ e^l - elf'^ ' 

(4) 



Fig. 3. — Migration history for a system that in the ab- 
sence of tides would have reached a minimum distance of 0.0004 
AU (0.1-Rq) and been declared lost; in the presences of tides it 
reaches a minimum distance of 0.013 AU and is later circularized at 
ap = 0.026 AU. The four panels are: top-left, planet eccentricity as 
a function of time (in years); top-right, relative inclination between 
the two orbit normals; bottom-left, planet semi-major axis (solid, 
in AU) and periapse (dotted); bottom-right, planet total orbital 
angular momentum (J, solid) and its component along the orbit 
normal of the binary (J^, dotted), both in arbitrary units. Kozai 
oscillation (which conserves Jz) has proceeded for barely half a 
cycle before the orbital energy of the planet is significantly dissi- 
pated and the planet is removed from the infiuence of the binary 
companion. Tidal dissipation operates afterwards (during which J 
is conserved). The inclination angle evolves little in this example. 

where ric = 27r/Pc, Up = 27r/Pp. The rate of tidal ec- 
centricity damping depends strongly on the periapse dis- 
tance. Consideri ng only tides raised on the planet, we 
obtain (HnHlMlD 



1 dCp 
e„ dt 



27kpGmp 1 
2RlQpnp q 



Rr 



1 



^pJ (l-e2)i3/2' 

where kp is the planet's tidal Love number, Qp its tidal 
dissipation factor and Rp its radius (see WM03). The 
mass ratio q = nip/m^,. Equating the two rates, we ob- 
tain the stalling periapse value. 



rp,st.n = 0.015AU 
'R. 



TO* 
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TO,, 
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Rj 



0.23M, 



TOc 



Qp 

\27QAUJ 



0.5 

1/6.5 



(6) 



where we have scaled variables by their representative 
values {Rj is the radius of Jupiter). Coincidentally, 
J^p, stall ^ tr, and it depends little on a variety of pa- 
rameters, including stellar mass, companion mass, planet 
mass, planet tidal Q factor, and planet initial orbit. This 
justifies our choices for these parameters in the numerical 
experiment. ^ In our simulation, most binaries that give 

^ This mechanism works for other types of planets like hot Nep- 
tuncs or super-earths. Substituting into Eq. |[B) values appropriate 
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Fig. 4. — Parameters for the binary systems that produce Kozai 
migration. In the top panel, the final ap (horizontal axis) is plotted 
against ac- Smaller values of the former are in general correlated 
with closer binaries (eq. |6]l, with most hot Jupiters arising from 
binaries with ac G [100, 1000] AU. The bottom left panel shows 
the distributions of initial (thin lines) and final (thicker lines) in- 
clinations between the two orbital planes - the solid curves include 
systems with ap < 0.025 AU, and the dotted curves all systems 
with final Op < 0.1 AU, (similarly in the right panel). The final 
inclination angles are much more sp read out, as the Kozai cycles 
convert inclination to eccentricity — IFabrvckv fc Tremaind ()2007l l 
gives a detailed explanation for the features. The bottom right 
panel shows the distribution of final angles ip between the stellar 
spin axis and the planetary orbit normal. Most systems (espe- 
cially the tightest ones) have ip < 50°, although some stars may 
spin retrograde relative to the planet orbit. 

rise to hot Jupiters have G [100, 1000] AU (Fig. H]) 
and we have scaled Oc here by roughly the median value. 
Tighter binaries are relatively unimportant - planets in 
many of these systems are dynamically unstable and are 
excluded from our study. 

In the subsequent tidal circularization, orbital angu- 
lar momentum is roughly conserved and the final ap ~ 
2rp,staii ~ 0.03 AU. 

The fraction of stars with Kozai migrated hot Jupiters 
is given by 

/<o.i ~ fb ■ fp ■ /k (7) 

where fb is the fraction of stars in binary systems, fp is 
the fraction of solar type stars with Jupiter mass planets 
formed at a few AU, and /Kozai is the fraction of plan- 
ets in binary star systems that un dergo Kozai migration 
to Op < 0.1 AU. Taking fh ^ 0.65 (iDuauennov fc Mayor! 
[Ml, fp ^ 0.07 (iMarcv et al.l l2005l). and f^.^.i ^ 0.025 
(this work), we suggest that, at a minimum, 10% of the 
known hot Jupiters may be due to Kozai migration. The 
most uncertain number is fp. The value of fp we have 
quoted is the observed fraction in the Keck sample, which 
is substantially complete up to Op « 3 AU. Assuming the 
number of planets per AU is flat up to ap = 30 AU gives 
fp = 0.12 and /<o.i = 0.002. There is some indication 
that the number of planets per AU is an increasing func- 
tor Neptune and Earth, we obtain similar values for stall- 



tion of flp. li fp — 0.5, more than half the hot Jupiters 
could be produced by the Kozai mechanism. 

3.2. Predictions of the Kozai Migration Scenario 

The number of hot Jupiters produced by Kozai mi- 
gration can be determined by observations in the near 
future, since Kozai migrated planets must have a num- 
ber of attributes. First, candidate Kozai hot Jupiters 
will reside in binary star systems, although the binary 
mass ratio may well be small; a brown dwarf companion 
can be dynamically as ef fective as a solar-type compan ion 
(eq. [6]). The study bv iDuquennov fc Mavoil (|1991h es- 
tablishes that ~ 60% of the stars in the solar neighbour- 
hood are actually binary or triple systems. While radial 
velocity surveys select against close binaries, studies by 
iRaghavan et al.l (|2006l ) show that at least 23% of radial 
velocity planet hosts have stellar companions. The dis- 
coveries of brown-dwarf companions to the planet bear- 
ing stars HP 3651 (M ugrauer et al. 2006) and HD 89744 
(jMugrauer et al.ll2004l) highlight the possibihty that the 
existence of dim companions will increase the known bi- 
nary fraction of planet bearing stars significantly. The 
Kozai scenario predicts that the binary fraction of hot 
Jupiters will be higher than that of systems with more 
distant planets. Binary-induced radial velocity trends in- 
duced on the primary by a stellar companion will be of 
order 

mc \ f 100 AUV 



5/ 



m/ s/yr, 



(8) 



where / is the sine of the angle between the line of 
sight and the stellar velocity. This is clearly detectable 
at the current sen sitivity of radial velocity surveys 
(|Wright et al.l[200l . The companion will also induce an 
astrometric acceleration of a few micro-arcsecond/yr/yr, 
detectable by SIM or GAIA. 

Second, Kozai systems have / G [30°, 150°], with / w 
90° not uncommon (Fig. 3]). In transiting systems the 
binary orbit will be in or near the plane of the sky. This 
can be tested via both radial velocity and astrometry. 

Third, the angle between the spin axis of the primary 
star (assumed to be the orbit normal of the planet at 
formation) and the present-day planet orbit normal will 
range from to 130 degrees (Fig. |4]) with the values 
between and 50 degrees being preferred. This angle 
can be determined if both the spin period of the star 
as well as its rotational velocity v sin i can be indepen- 
dently measured. The angle projected onto the plane of 
the sky, measurable using the Rossiter-McLaughlin ef- 
fect, will have a similar range. 

Fourth, the semimajor axis ratio Op/a'p with any second 
planet will be small. This results from the requirement 
that the precession rate ind uced by the second p lanet not 
break the Kozai resonance (iWu fc Murravl[200l . Radial 
velocity measurements can constrain the mass and semi- 
major axis of any nearby pl anetary companions to hot 
Jupiters (Wright et al. 2007). A corollary is that the 
fraction of multiplanet systems having hot Jupiters will 
be smaller than the fraction of single planet systems with 
hot Jupiters. 

Kozai-migrated planets dissipate many times their 
own binding energy during tidal circularization. 
lOgilvie fc LinI (|2004i) find that tidally dissipated energy 
is deposited throughout the bulk of the planet, raising 
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the possibility that the plane t will expand catastrophi- 
cally. In contrast, IWd ()2005[ ) concludes that energy is 
deposited exclusively near the photosphere, which would 
leave the planet intact. 

The theoretical situation is unclear, but the existence 
of hot Jupiters suggests an answer. A plot of Cp versus Op 
strongly suggests that the low Cp's of the hot Jupiters are 
the result of tidal circularization, as the observed Cp's fol- 
low closely the upper-bo und set by the tidal process (see, 
e.g.. Fig. 1 of 'Wu' '2003). If so, most or all hot Jupiters 
have experienced rapid tidal heating and survived. 

Another concern with the Kozai picture is raised by 
the Rossiter-McLaughlin measurement of stellar obliq- 
uity, cu rrently available for 5 transit ing planets (see Ta- 
ble 2 of lFabrvckv fc Tremainj|2007f) . All are consistent 
with zero obliquity. Taken at face value, this is at vari- 
ance with the above Kozai prediction. * 

3.3. Alternatives to Kozai Migration 

In Kozai migration, it is important that rp evolves on 
a time scale no shorter than the tidal precession time 
scale; if rp were to suddenly plunge from above to below 
the Roche radius, as for example would be the case if 
two planets suffered a close encounter, the inward scat- 
tered planet would not be stalled outside r^. Instead it 
would suffer rapid mass loss and likely be lost. In that 
case there will be a cut -off in the distribution of ap at 
2rB (|Ford fc Rasioll2OO60 . but not a pile-up. 



Migration in a gas disk may also produce hot Jupiters. 
If the disk extends all the way to the star, one would 
observe a cut-off at Op ~ rp- if the disk is truncated, 
e.g., by stellar magnetic fields ()Lin et al.lll996f) . a feature 
will appear at an orbital period half that of the inner 
edge of the disk. However, spin periods and magnetic 
fields of accreting stars show a substantial dispersion, 
which would lead to a rather broad distribution in the 
disk inner radii, and hence a smeared out feature in the 
distribution of planetary semimajor axis. 

We have studied the role of a binary companion in 
increasing Cp and causing a gradual collapse in rp. 
But it is also plausible that so ft planet-planet scatter- 
ing can gradually dec rease rp (jJuric fc Tremaind l2007t 
IChatteriee et al.l(2007f) . Moreover, Kozai oscillations can 
also be excited by a second planet,^ in the absence of a 
binary stellar companion. As long as these or other pro- 
cesses produce gentle eccentricity driving on 10^ to 10* 
year timescales, tidal effects will halt the periapse evolu- 
tion when rp ^ r^. Tidal circularization then pushes the 
planets out to Op ~ 2rfi and produces a narrow pile-up 
of hot Jupiters there. 



We acknowledge helpful discussions with Scott Gaudi, 
Daniel Fabrycky and Andrew Gould, as well as NSERC 
discovery grants to YW and NM, and an NSERC under- 
graduate fellowship to MR (summer 2006). 



* HD147506 l lWinii et al.l l2007i : [Loeillct ct al. 2007'), a 1.3Mq 
star with a massive planet, may have experienced tidal synchro- 
nization in its surface-layer that would alter its apparent rotation 



axis 
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This second planet can be placed on a highly inclined orbit by, 



e.g., planet-planet scattering. 
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